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Abstract 

We derive the integral fluctuation theorem around a nonequilibrium stationary state for friction- 
less and soft core granular particles under an external vibration achieved by a balance between 
an external vibration and inelastic collisions. We also discuss the connection between the integral 
fluctuation theorem and the generalized Green-Kubo formula. 
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I. INTRODUCTION 



One of the most remarkable achievements in recent nonequihbrium statitistical mechanics 
is to demonstrate the existence of some nonequihbrium identities such as the generahzed 
Green-Kubo relation [l, 2|, various fluctuation theorems js Sl and the Jarzynski equality 
as well as the mutual relationship 10| . These identities are exact and reproduce the con- 
ventional Green-Kubo formula, the second law of thermodynamics and Onsager's reciprocal 
relation in specific limits. Therefore, these identities are regarded as fundamental relations 
in nonequihbrium statistical mechanics. 

Although it has been believed that these identities are supported by the local time- 
reversal symmetry or the detailed balance condition, some experiments suggest the exis- 
tence of fluctuation theorem or related equation even in granular systems which do not 

ll-|l6|, t 



have any time reversal symmetry 



lough there exists a counter argument 



remarkable that Puglisi and his coworkers 



17|. It is 



18l-l2l| clarified that granular fluids do not hold 



the conventional fluctuation theorem but have only the second type fluctuation theorem 



by Evans and Searles 



221 ]. As long as the author's knowledge, however, there is only a 



paper by Chong et al. which has proven the existence of both the generalized Green-Kubo 
relation and the integral fluctuation theorem 8|] for a granular system under a steady plane 



shear 



12| . They also developed the representation of a nonequilbrium steady-state distribu- 



tion function 



23| and the liquid theory for sheared dense granular systems. 2J]. Recently, 



the present author 



25[ extended their previous formulation to discuss the nonlinear response 



theory around a nonequihbrium steady state, and also demonstrate that the conventional 
fluctuation theorem can be obtained as well as the Jarzynski equality. 



12 



are only valid for sheared cases. 



11 



13 
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In this paper. 



Unfortunately, our previous theoretical studies 
but the most of experiments adopt vibrating granular gases 
thus, we derive the Jarzynski equality for soft core granular gases under vibrations around 
a nonequihbrium steady state. We also discuss the relationship between the gerenalized 
Green-Kubo formula and the integral fluctuation theorem. 

The organization of this paper is as follows. In Sec. II, we summarize the general 
framework of Liouville equation and some identities which are used in this paper. Section 
III which is the main part of this paper consists of two parts. In the first part (Sec. Ill A) 
we discuss the the derivation of the integral fluctuation theorem (IFT). In the second part 



2 



(Sec. Ill B) we also derive a standard fluctuation theorem and the generalized Green-Kubo 
formula from IFT. In Sec. IV we discuss our results and we give conclusion in section V. In 
Appendix A, we briefly summarize some operators' identities. 



II. LIOUVILLE EQUATION 

Let us consider a system of identical soft spherical and smooth dissipative particles. 
We assume that particles are mono dispersed, which are characterized by their diameter d 
and the mass m. The particles are influenced by the gravity with the acceleration constant g 
in z— direction If we use a box fixed frame, each particle feels the acceleration —g+Auj'^ cos ut 
in 2;— direction with the amplitude A and the angular acceleration u. Moreover, we should 
introduce a confined potential that particles cannot penetrate the bottom plate. 

The basic equation for the statistical mechanics of frictionless granular particles under 

nnn 

such a vibrations is the Liouville equation. |26H32| The argument in this section is parallel 
to that in Ref.[l|. Let iC{t) be the total Liouvillian which operates an arbitrary function 
A{T{t)) starting from t = as 

dA{T{t)) 
dt 

where 



f/^(0, t)iC{t)A{T), A{T{t)) = U^iO, t)A{r), (1) 



= I '^^^ I ds2--- dSniC{sn) ■ ■ ■iC{s2)iC{si), (2) 

, n Jo Jo Jo 



n=0 ■ 

iN 



and T{t) = {ri{t),Pi{t)}iLi with the abbreviation F = r(0). We note that there are some 
trivial relations for U^{to,t) such as 

U^{to,t) = U^ito,s)U^{s,ty, U^ito,t)fiT{to)) = Km) (3) 

for an arbitrary function f(T(t)). 

The total Liouvillian consists of three parts, the elastic part, the viscous part and the 
part from an external vibration. We write iC{t) as 

iC{t) = zi:(^')(r) + z£("^^)(r) + z£^"*(r,t), (4) 
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where iC^^^^ is the elastic colhsion part, 



N 



1=1 



m dvi ' dp 



Here, we assume that the elastic force can be represented by the summation of the pairwise 
force f/^') = Fjf with 

= -^-^ = Qid - n,)f{d - r,,)r,„ (6) 

where we have introduced the pair-wise potential u{rij), = fi—fj, fij = |rjj|, fij = Vij/vij, 
and the Heviside function 0(x) which satisfies B(x) = 1 for x > and 0(x) = for 
otherwise. The elastic repulsive force f{x) is proportional to x for the linear spring model, 
and to x"^/^ for the Hertzian contact model. 

Similarly, the viscous Liouvillian is the contribution of inelastic collisions: 

N 

p(vis) _ 

dpi 



,£(vis)(r)^^^(vis).^^ (7) 



i=l 

where F^^^^^ is the viscous force acting on i—th particle represented by i^T*-^"^-* = ^^^^F^J^^^ 
with 

F.^™^ = -rijQ{d - r,j)C{d - rij){Vij ■ Vij). 

= -fijj^{rij){vij - fij). (8) 

Here we have introduced Vij = Vij = dvi/dt, and 

J^{r) = e{d - r)C{d - r) (9) 

with the viscous function which is typically a constant or ({x) oc x^^^. The Liouville 
operator representing the vibration iC^^^{t) is given by 

N 

dpi 

where the vibrating force is given by 



iC^^\T,t) = J2Ft'\t).^, (10) 



f/^-*) {t) = z\ m{-g + Au^ cos ut) - ^^^f^ \ = zFt""'^ (t) (11) 
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in a box fixed frame, wliere z is tfie unit vector in z direction, and Vcyxi^z) represents a 
confined potential in a box such as 



Kxt(2) = Kjexp {-zjQ^ 



(12) 



to prevent grains from penerating tlie bottom plate of tlie container. 

It should be noted that the Liouvillian is not self-adjoint, because of the violation of 
time-reversal symmetry for each collision. The adjoint Liouvillian is defined through the 
equation of the phase function or the A^— body distribution function p(r, t) 



p(r, t) = u^{t, o)p(r, 0), = -^c\t)p{T, t), 

where 

.t 

n=0 

The adjoint Liouvillian satisfies 

tC\r,t) = iC(T,t)+A{T) 



(13) 



t rsi 

1 / dS2 
^0 



(14) 



where 



A(r) = ^ •r(r) 



(15) 



(16) 



is the phase volume contraction. We note that A(r) in our system does not depend on t 
explicitly, which can be written as 



I]) 



(17) 



for t > 0. The phase volume contraction A(r) is directly related to the change of Jacobian 

dT{t) r rt 



dT 



exp 



drA(r(r)) 



(18) 



where A(r(t)) = f/-^(0, t)A(r)f/^(t, 0). Note that the time evolution of an arbitrary physical 
function A{V{t)) is given by A{V{t)) = U^{0,t)A(T)U^{t,0), where we have introduced 
U^{t, 0) = exp[— i Jq* dsC{s)]. 

The average of a physical quantity is defined as 



{A{rm ^ I dTp{T,o)A{T{t)) = I drA{r)p{T,t). 



(19) 



From Eqs. ([T]), (fT3|) and (fT9|) we obtain the relations 

J dTp{T)U^{0,t)A{T) = J dTA{T)U^{t,0)p{T,0) (20) 

and 

j dTp{T)iC{t)A{T) = - j dTA{T)ia{t)p{T,Q). (21) 

In the last part of this section, let us introduce a stationary distribution to characterize the 
quasi-periodic motion of granular particles under the periodic vibration. In the stationary 
process, the initial distribution function p(r, 0) may have the form 

e-^o(r) 

p(r,o) = pini(r) = ^^, (22) 

where /o(r) = /(F, t = 2mr/uj) with an arbitrary integer n. and Z = j dTe~^°^^\ Note that 
/o(r) is an arbitrary function of F, and thus, this choice is quite general for the argument. In 
the stationary process, we assume that an average of an arbitray function A{T{t)) satisfies 
the perodic condition 

A{T(t+'^-^\t\ = {A{V{t)) (23) 



for any nonzero integer n. This assumption is reasonable because of the periodicity of the 
Liouvillian: 

iC (t + = iC{t) (24) 



for an arbitrary integer n. 



III. FLUCTUATION THEOREM AND GREEN-KUBO FORMULA 

This section consists of two parts. The first part is dedicated to the derivation of the 
integral fluctuation theorem (IFT). The second part discusses the derivations of both the 
standard fluctuation theorem and the generalized Green-Kubo formula. 



A. Integral fluctuation theorem 



The integral fluctuation theorem (IFT) is a kind of the fluctuation theorem, which is 
related to Jarzynski equality, jsl The relation between Jarzynski equality and the fluctuation 



theorem has been investigated extensively. 



10| Although IFT is an important identity for 
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the stationary sheared granular systems [l2|], and this equahty plays a fundamental role even 
in granular systems under a vibration. 

To demonstrate the existence of the IFT, we consider a system characterized by the 
following time-dependent Hamilitonian : 



^o(r(t))=$:^+^5:.Mt)). 



(25) 



We also assume that the initial condition satisfies the canonical distribution 

e-/3^o(r) 



Peq(r) 



(26) 



where /3 is the inverse temperature and Z{/3) = / dTe 

In this case, it is easy to verify the conservation of the normalization factor, i.e. Z{I3) 
J rfre-^^o(^) = / rfr(t)e-'^^o(^W. Because the time derivative of Ho{T{t)) is given by 

Mm) = J2''^At)Ft'\t) - 27^(^(t)) 

i 

with Hq = dHo/dt, Vi{t) = dri{t)/dt, 



(27) 



T{vis) 



2 

I] J ! 



(28) 



and Vij = Vi — Vj, the conservation of Z{I3) satisfies 

g-/3Ho(r(t)) 



dT{t)- 



dr 



dm 



dT 



g-/3Ho(r) 



exp 



/3 / dr{J2^^.{r)Ft'\r)-2nmm 



g-/3//o(r) 
dT — TTT^. — exp 



ZW) 



rfrfieq(r(r)) 



exp 



Meq(r(r)) 



eq 



where we have used Eq. (fT8!) for the third equality, and introduced 

neq(r(t)) ^ /3 5^i;,.(t)Ff^*^(t) - 2/37^(^(t)) - A(r(t)). 



(29) 



(30) 



Equation fl2U]) associated with Eq. fl5Ul) is the IFT for granular fluids under the vibration. 

It is a characteristic feature for dissipative systems that the phase volume contraction 
A(r) is involved in r2cq(r(t)). Thus, the right hand side of the Jarzynski equality (!29|) for 
dissipative cases cannot be represented by the work done by the external force. 



The IFT (129|) is directly reduced to an inequality 



rfr(l]eq(r(r))) > 



(31) 



with the aid of Jenssen's inequality. 

Now, let us extend the IFT to the case of starting from an arbitrary distribution pini(r). 
In this case, the IFT can be rewritten as 

p-io(m) 



1 



dT{t)- 



Z 



dT 



dm 



dT 



-Mr) 



Z 



exp 



rfr/o(r(r)) 



e--fo{r) 
dr — - — exp 



Z 



drfi(r(r)) 



exp 



rfrl](r(r)) 



(32) 



where fi(r(t)) = Jo(r(t)) - A(r(t)), and Z = / dFe-^o^^) = / rfr(t)e-^«(rW). Equation ([32 
is also reduced to 



cir(fi(r(r),r)) > 0. 



(33) 



B. The standard Fluctuation Theorem and GerenaHzed Green-Kubo formula 



The direct consequences of Eq. fl32|) are two important relations, the conventional fluctu- 
ation theorem and the generalized Green-Kubo formula, from Eq. (l32l) . Here, let us illustrate 
how to derive these two relations. 

It is straightforward to derive the standard fluctuation theorem from IFT [l| , if we assume 
the time reversal symmetry in pini(r) as 



Pini(r^) = pi,i(r) 



(34) 



where represents the time reversal operation of F, i.e. F^ = {r^, — pj}^]^. Let us derive 
the standard fluctuation theorem from the integral fluctuation theorem ( 132|) . Equation ( 132|1 
can be rewritten as 

'rfFpM(r)e-*^ = y rfF^pM(r^), (35) 

where we have introduced fi^ = j (irr2(F(r)) for t = 2mi/uj for an arbitrary positive inte- 
ger n and used the fact that the normalization is unchanged for J dT'^ pi^i{T^) = 1. With the 



aid of Eq.dalD and dT'^ = dT\dr^/dT\, Eq. ([35]) means that the Jacobian {dV^/dri is given 
by exp[—t^lt]. Therefore, we can write the probabihty of fit = A for fi^ = ^ dr^l{T(r)) 
and n{T) = fi(r^) as 



Prob(fii = A) 



dT^pUT^)5mT^)-A) 
drp,^i{T)e~'^'6iWt{T) - A) 



-At 



e-^*Prob(a = A), 



rfrpi„i(r)5(fit(r) - A) 



(36) 



where we have used fi(r"^) = Q{{V^}'^) = Q(T) in the second hne. Note that our fluctuation 
theorem 

Prob(a = A) = e-^*Prob(n; = A) (37) 

differs from the standard form, because our system does not have the time reversal symmetry. 

Next, let us derive the generalized Green-Kubo formula from Eq. fl32l) following the argu- 
ment in Ref.(l2|. Equation leads to the relation rfrpini(r)e-^o ^^^^(rM) = rfr(t)pini(r(t)) 
in a co-moving frame. Multiplying e-^o dTn{T{T)) ^^j. side, we obtain 

rfrpi„i(r) = c/r(t)pini(r(t))exp 

r 

drpini(r) exp 
rffpini(f ) exp 



rfrfi(r(r)) 
L"'o 

rfrfi(f(r)) 
rfrfi(f(-r)) 



(38) 



where we have introduced T = r{t) and used r(r) = T(t + t). Therefore, we reach the 
conservation of the probability: 

(e^o-^-f^Cn-))) = 1. (39) 



From Eq.f p9|) . we obtain the time evoluation of the distribution function: 

p(r,t) = e/o'^-^(r(-))p^^^(r). 



(40) 



The generalized Green-Kubo formula is the direct consequence of the IFT. The differen- 
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tiation of Eq. (fT9|) with the help of Eq. (HO|) . we obtain 
d 



^^{A{Tm = J dTAiTMTi-t))piT,t) = J dTU^{t,0){A{T{t)MT)}piT,t) 

= J dTA{T{t))n{T)U^{0,t)p{T,t) = J dTA{T{t))n{T)e^o'i-Mn-))p(T{t),t) 
= J dTA{T{t))n{T)^^ = {A{T{tmT)). (41) 
This equation can be integrated over t as 

(A(r(t))) = {A{T)) + fds{A{T{smT)). (42) 
This result depends on pini(r). Therefore, the formal response theory can be written as 

POO 

(M(r)) = / dtiAiTitmr)), (43) 

Jo 

where 6A{T) = linit^oo ^(r(t)) - A{T). 
IV. DISCUSSION 

In this paper, we obtain exact nonequilibrium relations. To verify the validity, we may 
need numerical simulations. It should be noted that the verification of the generalized 
Green-Kubo formula is not difficult by the direct simulation, but the confirmation of the 
integral fluctuation theorem by simulations is not easy because of the limitation of numerical 

ization condition fj39|) could not be achieved in the simulation 



accuracy. Indeed, the norma 



even for thermostat systems. [34] We should stress that the generalized Green-Kubo formula 
can be used for dense granular systems above the jamming transition. 

Evan when we confirm the validity of the nonequilibrium relations such as the generalized 
Green-Kubo formula, it is not easy to calculate the correlation function Eq. (l42|) or Eq.( l43|l . 
One of possible methods is to use the mode-coupling theory (MOT). It is helpful to a pply 



29 



36 



38| 



MCT for granular liquids to characterize theology near the jamming transition. 
It is notable that Ref. 3JJ] develops a linear response theory for a sheared thermostat system 
around a nonequilibrium steady state. The application of this method will be discussed 
elsewhere. 



10 



V. SUMMARY 



We have developed some exact relations for frictionless granular fluids under vibrations. 
We derived the integral fluctuation theorem around a nonequilibrium steady state. We 
also derived the formal expression of the distribution function. We flnally obtained the 
generalized Green-Kubo formula around a nonequilibrium steady state. 

In this paper, we focus on the detailed analytic calculation on the granular fluids under 
the vibration. The systematic check in terms of the simulations will be reported elsewhere. 
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Appendix A: Some operators' identities 

Let us consider the time evolution of r(t; to) deflned by 

r(t;to) = f/-^(to,t)r(to), (Al) 

where we have explicitly written the initial time to. 

By using U^{to,t) and U^{t,to) it is notable that there is an important relation for 
U^{to,t): 

AiTit)) = U^{to,t)AiT{to))U^{t,to) = U^{to,t)A{r{to)). (A2) 
The proof of (lA2p is straightforward. The right hand side of Eq. (IA2p can be rewritten as 

U^ito,t)A{T{to)) = U^ito,t)A{T{to))U^{t,to)U^it,to)l = U^{to,t)A{T{to))U^it,to)l 
= U^{to,t)A{T{to))U^{t,to), (A3) 

where we have used U^{to,t)l = 1 for a constatnt 1. When we use Eq. ( 1A2I) . we readily 
obtain 

U^{to, t)A{T{to))B{T{to)) = A{T{t- to)) ■ B{T{t; to)). (A4) 
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Indeed, the left hand side of this equation can be rewritten as 

f/^(to, t)A{T{to))B{r{to)) = U^ito, t)A{r{to))UUt, to)U^ito, t)5(r(to))f/^(t, to) 

= A{T{t;to))-B{T{t;to)), (A5) 

which is the end of the proof of Eq. ( 1A4|) . 
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